Evolution of a Human-Specific De Novo Open Reading Frame and Its Linked Transcriptional Silencer.

In the human genome, two short open reading frames (ORFs) separated by a transcriptional silencer and a small intervening sequence stem from the gene SMIM45. The two ORFs show different translational characteristics, and they also show divergent patterns of evolutionary development. The studies presented here describe the evolution of the components of SMIM45. One ORF consists of an ultra-conserved 68 amino acid (aa) sequence, whose origins can be traced beyond the evolutionary age of divergence of the elephant shark, ~462 MYA. The silencer also has ancient origins, but it has a complex and divergent pattern of evolutionary formation, as it overlaps both at the 68 aa ORF and the intervening sequence. The other ORF consists of 107 aa. It develops during primate evolution but is found to originate de novo from an ancestral non-coding genomic region with root origins within the Afrothere clade of placental mammals, whose evolutionary age of divergence is ~99 MYA. The formation of the complete 107 aa ORF during primate evolution is outlined, whereby sequence development is found to occur through biased mutations, with disruptive random mutations that also occur but lead to a dead-end. The 107 aa ORF is of particular significance, as there is evidence to suggest it is a protein that may function in human brain development. Its evolutionary formation presents a view of a human-specific ORF and its linked silencer that were predetermined in non-primate ancestral species. The genomic position of the silencer offers interesting possibilities for the regulation of transcription of the 107 aa ORF. A hypothesis is presented with respect to possible spatiotemporal expression of the 107 aa ORF in embryonic tissues.

ProTInSeq: transposon insertion tracking by ultra-deep DNA sequencing to identify translated large and small ORFs.

Identifying open reading frames (ORFs) being translated is not a trivial task. ProTInSeq is a technique designed to characterize proteomes by sequencing transposon insertions engineered to express a selection marker when they occur in-frame within a protein-coding gene. In the bacterium Mycoplasma pneumoniae, ProTInSeq identifies 83% of its annotated proteins, along with 5 proteins and 153 small ORF-encoded proteins (SEPs; ≤100 aa) that were not previously annotated. Moreover, ProTInSeq can be utilized for detecting translational noise, as well as for relative quantification and transmembrane topology estimation of fitness and non-essential proteins. By integrating various identification approaches, the number of initially annotated SEPs in this bacterium increases from 27 to 329, with a quarter of them predicted to possess antimicrobial potential. Herein, we describe a methodology complementary to Ribo-Seq and mass spectroscopy that can identify SEPs while providing other insights in a proteome with a flexible and cost-effective DNA ultra-deep sequencing approach.

Blastocrithidia nonstop mitochondrial genome and its expression are remarkably insulated from nuclear codon reassignment.

The canonical stop codons of the nuclear genome of the trypanosomatid Blastocrithidia nonstop are recoded. Here, we investigated the effect of this recoding on the mitochondrial genome and gene expression. Trypanosomatids possess a single mitochondrion and protein-coding transcripts of this genome require RNA editing in order to generate open reading frames of many transcripts encoded as 'cryptogenes'. Small RNAs that can number in the hundreds direct editing and produce a mitochondrial transcriptome of unusual complexity. We find B. nonstop to have a typical trypanosomatid mitochondrial genetic code, which presumably requires the mitochondrion to disable utilization of the two nucleus-encoded suppressor tRNAs, which appear to be imported into the organelle. Alterations of the protein factors responsible for mRNA editing were also documented, but they have likely originated from sources other than B. nonstop nuclear genome recoding. The population of guide RNAs directing editing is minimal, yet virtually all genes for the plethora of known editing factors are still present. Most intriguingly, despite lacking complex I cryptogene guide RNAs, these cryptogene transcripts are stochastically edited to high levels.

cirCodAn: A GHMM-based tool for accurate prediction of coding regions in circRNA.

Studies focusing on characterizing circRNAs with the potential to translate into peptides are quickly advancing. It is helping to elucidate the roles played by circRNAs in several biological processes, especially in the emergence and development of diseases. While various tools are accessible for predicting coding regions within linear sequences, none have demonstrated accurate open reading frame detection in circular sequences, such as circRNAs. Here, we present cirCodAn, a novel tool designed to predict coding regions in circRNAs. We evaluated the performance of cirCodAn using datasets of circRNAs with strong translation evidence and showed that cirCodAn outperformed the other tools available to perform a similar task. Our findings demonstrate the applicability of cirCodAn to identify coding regions in circRNAs, which reveals the potential of use of cirCodAn in future research focusing on elucidating the biological roles of circRNAs and their encoded proteins. cirCodAn is freely available at https://github.com/denilsonfbar/cirCodAn.

Widespread stable noncanonical peptides identified by integrated analyses of ribosome profiling and ORF features.

Studies have revealed dozens of functional peptides in putative 'noncoding' regions and raised the question of how many proteins are encoded by noncanonical open reading frames (ORFs). Here, we comprehensively annotate genome-wide translated ORFs across five eukaryotes (human, mouse, zebrafish, worm, and yeast) by analyzing ribosome profiling data. We develop a logistic regression model named PepScore based on ORF features (expected length, encoded domain, and conservation) to calculate the probability that the encoded peptide is stable in humans. Systematic ectopic expression validates PepScore and shows that stable complex-associating microproteins can be encoded in 5'/3' untranslated regions and overlapping coding regions of mRNAs besides annotated noncoding RNAs. Stable noncanonical proteins follow conventional rules and localize to different subcellular compartments. Inhibition of proteasomal/lysosomal degradation pathways can stabilize some peptides especially those with moderate PepScores, but cannot rescue the expression of short ones with low PepScores suggesting they are directly degraded by cellular proteases. The majority of human noncanonical peptides with high PepScores show longer lengths but low conservation across species/mammals, and hundreds contain trait-associated genetic variants. Our study presents a statistical framework to identify stable noncanonical peptides in the genome and provides a valuable resource for functional characterization of noncanonical translation during development and disease.

Deep whole-genome analysis of 494 hepatocellular carcinomas.

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

Approaches to probe and perturb long noncoding RNA functions in diseases.

Long noncoding RNAs (lncRNAs) are a class of RNA molecules exceeding 200 nucleotides in length that lack long open-reading frames. Transcribed predominantly by RNA polymerase II (>500nt), lncRNAs can undergo splicing and are produced from various regions of the genome, including intergenic regions, introns, and in antisense orientation to protein-coding genes. Aberrations in lncRNA expression or function have been associated with a wide variety of diseases, including cancer, cardiovascular diseases, diabetes, and neurodegeneration. Despite the growing recognition of select lncRNAs as key players in cellular processes and diseases, several challenges obscure a comprehensive understanding of their functional landscape. Recent technological innovations, such as in sequencing, affinity-based techniques, imaging, and RNA perturbation, have advanced functional characterization and mechanistic understanding of disease-associated lncRNAs.

Sequences at gene segment termini inclusive of untranslated regions and partial open reading frames play a critical role in mammalian orthoreovirus S gene packaging.

Mammalian orthoreovirus (MRV) is a prototypic member of the Spinareoviridae family and has ten double-stranded RNA segments. One copy of each segment must be faithfully packaged into the mature virion, and prior literature suggests that nucleotides (nts) at the terminal ends of each gene likely facilitate their packaging. However, little is known about the precise packaging sequences required or how the packaging process is coordinated. Using a novel approach, we have determined that 200 nts at each terminus, inclusive of untranslated regions (UTR) and parts of the open reading frame (ORF), are sufficient for packaging S gene segments (S1-S4) individually and together into replicating virus. Further, we mapped the minimal sequences required for packaging the S1 gene segment into a replicating virus to 25 5' nts and 50 3' nts. The S1 UTRs, while not sufficient, were necessary for efficient packaging, as mutations of the 5' or 3' UTRs led to a complete loss of virus recovery. Using a second novel assay, we determined that 50 5' nts and 50 3' nts of S1 are sufficient to package a non-viral gene segment into MRV. The 5' and 3' termini of the S1 gene are predicted to form a panhandle structure and specific mutations within the stem of the predicted panhandle region led to a significant decrease in viral recovery. Additionally, mutation of six nts that are conserved across the three major serotypes of MRV that are predicted to form an unpaired loop in the S1 3' UTR, led to a complete loss of viral recovery. Overall, our data provide strong experimental proof that MRV packaging signals lie at the terminal ends of the S gene segments and offer support that the sequence requirements for efficient packaging of the S1 segment include a predicted panhandle structure and specific sequences within an unpaired loop in the 3' UTR.

Translation of non-canonical open reading frames as a cancer cell survival mechanism in childhood medulloblastoma.

A hallmark of high-risk childhood medulloblastoma is the dysregulation of RNA translation. Currently, it is unknown whether medulloblastoma dysregulates the translation of putatively oncogenic non-canonical open reading frames (ORFs). To address this question, we performed ribosome profiling of 32 medulloblastoma tissues and cell lines and observed widespread non-canonical ORF translation. We then developed a stepwise approach using multiple CRISPR-Cas9 screens to elucidate non-canonical ORFs and putative microproteins implicated in medulloblastoma cell survival. We determined that multiple lncRNA-ORFs and upstream ORFs (uORFs) exhibited selective functionality independent of main coding sequences. A microprotein encoded by one of these ORFs, ASNSD1-uORF or ASDURF, was upregulated, associated with MYC-family oncogenes, and promoted medulloblastoma cell survival through engagement with the prefoldin-like chaperone complex. Our findings underscore the fundamental importance of non-canonical ORF translation in medulloblastoma and provide a rationale to include these ORFs in future studies seeking to define new cancer targets.

Quantitative comparison of nuclear transport inhibition by SARS coronavirus ORF6 reveals the importance of oligomerization.

Open Reading Frame 6 (ORF6) proteins, which are unique to severe acute respiratory syndrome-related (SARS) coronavirus, inhibit the classical nuclear import pathway to antagonize host antiviral responses. Several alternative models were proposed to explain the inhibitory function of ORF6 [H. Xia et al., Cell Rep. 33, 108234 (2020); L. Miorin et al., Proc. Natl. Acad. Sci. U.S.A. 117, 28344-28354 (2020); and M. Frieman et al., J. Virol. 81, 9812-9824 (2007)]. To distinguish these models and build quantitative understanding of ORF6 function, we developed a method for scoring both ORF6 concentration and functional effect in single living cells. We combined quantification of untagged ORF6 expression level in single cells with optogenetics-based measurement of nuclear transport kinetics, using methods that could be adapted to measure concentration-dependent effects of any untagged protein. We found that SARS-CoV-2 ORF6 is ~15 times more potent than SARS-CoV-1 ORF6 in inhibiting nuclear import and export, due to differences in the C-terminal region that is required for the NUP98-RAE1 binding. The N-terminal region was required for transport inhibition. This region binds membranes but could be replaced by synthetic constructs which forced oligomerization in solution, suggesting its primary function is oligomerization. We propose that the hydrophobic N-terminal region drives oligomerization of ORF6 to multivalently cross-link the NUP98-RAE1 complexes at the nuclear pore complex, and this multivalent binding inhibits bidirectional transport.

RNA-binding protein Nocte regulates Drosophila development by promoting translation reinitiation on mRNAs with long upstream open reading frames.

RNA-binding proteins (RBPs) with intrinsically disordered regions (IDRs) are linked to multiple human disorders, but their mechanisms of action remain unclear. Here, we report that one such protein, Nocte, is essential for Drosophila eye development by regulating a critical gene expression cascade at translational level. Knockout of nocte in flies leads to lethality, and its eye-specific depletion impairs eye size and morphology. Nocte preferentially enhances translation of mRNAs with long upstream open reading frames (uORFs). One of the key Nocte targets, glass mRNA, encodes a transcription factor critical for differentiation of photoreceptor neurons and accessory cells, and re-expression of Glass largely rescued the eye defects caused by Nocte depletion. Mechanistically, Nocte counteracts long uORF-mediated translational suppression by promoting translation reinitiation downstream of the uORF. Nocte interacts with translation factors eIF3 and Rack1 through its BAT2 domain, and a Nocte mutant lacking this domain fails to promote translation of glass mRNA. Notably, de novo mutations of human orthologs of Nocte have been detected in schizophrenia patients. Our data suggest that Nocte family of proteins can promote translation reinitiation to overcome long uORFs-mediated translational suppression, and disruption of this function can lead to developmental defects and neurological disorders.

Improved super-resolution ribosome profiling reveals prevalent translation of upstream ORFs and small ORFs in Arabidopsis.

A crucial step in functional genomics is identifying actively translated ORFs and linking them to biological functions. The challenge lies in identifying short ORFs, as their identification is greatly influenced by data quality and depth. Here, we improved the coverage of super-resolution Ribo-seq in Arabidopsis (Arabidopsis thaliana), revealing uncharacterized translation events for nuclear, chloroplastic, and mitochondrial genes. Assisted by a transcriptome assembly, we identified 7,751 unconventional translation events, comprising 6,996 upstream ORFs (uORFs) and 209 downstream ORFs on annotated protein-coding genes, as well as 546 ORFs in presumed noncoding RNAs. Proteomic data confirmed the production of stable proteins from some of these unannotated translation events. We present evidence of active translation from primary transcripts of trans-acting small interfering RNAs (TAS1-4) and microRNAs (pri-MIR163 and pri-MIR169) and periodic ribosome stalling supporting cotranslational decay. Additionally, we developed a method for identifying extremely short uORFs, including 370 minimum uORFs (AUG-stop), and 2,921 tiny uORFs (2 to 10 amino acids) and 681 uORFs that overlap with each other. Remarkably, these short uORFs exhibit strong translational repression as do longer uORFs. We also systematically discovered 594 uORFs regulated by alternative splicing, suggesting widespread isoform-specific translational control. Finally, these prevalent uORFs are associated with numerous important pathways. In summary, our improved Arabidopsis translational landscape provides valuable resources to study gene expression regulation.

A method for evaluating of RNA's coding potential using the interaction effects of open reading frames and high-energy scalograms.

The identification and function determination of long non-coding RNAs (lncRNAs) can help to better understand the transcriptional regulation in both normal development and disease pathology, thereby demanding methods to distinguish them from protein-coding (pcRNAs) after obtaining sequencing data. Many algorithms based on the statistical, structural, physical, and chemical properties of the sequences have been developed for evaluating the coding potential of RNA to distinguish them. In order to design common features that do not rely on hyperparameter tuning and optimization and are evaluated accurately, we designed a series of features from the effects of open reading frames (ORFs) on their mutual interactions and with the electrical intensity of sequence sites to further improve the screening accuracy. Finally, the single model constructed from our designed features meets the strong classifier criteria, where the accuracy is between 82% and 89%, and the prediction accuracy of the model constructed after combining the auxiliary features equal to or exceed some best classification tools. Moreover, our method does not require special hyper-parameter tuning operations and is species insensitive compared to other methods, which means this method can be easily applied to a wide range of species. Also, we find some correlations between the features, which provides some reference for follow-up studies.

Small protein plays with big networks.

Thousands of small proteins, called microproteins, are encoded in small open reading frames (smORFs) throughout the genome. Despite assumptions that these proteins would be too small to properly fold and function, a recent study by Chen et al. identifies the surprisingly complex roles of one such microprotein.

The interplay between cis- and trans-acting factors drives selective mRNA translation initiation in eukaryotes.

Translation initiation consists in the assembly of the small and large ribosomal subunits on the start codon. This important step directly modulates the general proteome in living cells. Recently, genome wide studies revealed unexpected translation initiation events from unsuspected novel open reading frames resulting in the synthesis of a so-called 'dark proteome'. Indeed, the identification of the start codon by the translation machinery is a critical step that defines the translational landscape of the cell. Therefore, translation initiation is a highly regulated process in all organisms. In this review, we focus on the various cis- and trans-acting factors that rule the regulation of translation initiation in eukaryotes. Recent discoveries have shown that the guidance of the translation machinery for the choice of the start codon require sophisticated molecular mechanisms. In particular, the 5'UTR and the coding sequences contain cis-acting elements that trigger the use of AUG codons but also non-AUG codons to initiate protein synthesis. The use of these alternative start codons is also largely influenced by numerous trans-acting elements that drive selective mRNA translation in response to environmental changes.

A small stretch of poor codon usage at the beginning of dengue virus open reading frame may act as a translational checkpoint.

OBJECTIVE: Rare codons were previously shown to be enriched at the beginning of the dengue virus (DENV) open reading frame. However, the role of rare codons in regulating translation efficiency and replication of DENV remains unclear. The present study aims to clarify the significance of rare codon usage at the beginning of DENV transcripts using the codon adaptation index (CAI). METHODOLOGY: CAIs of the whole starting regions of DENV transcripts as well as 18-codon sliding windows of the regions were analyzed. RESULTS: One of the intriguing findings is that those rare codons do not typically result in uniformly low CAI in the starting region with rare codons. However, it shows a notable local drop in CAI around the 50th codon in all dengue serotypes. This suggests that there may be a translational checkpoint at this site and that the rare codon usage upstream to this checkpoint may not be related to translational control.

Ryegrass mottle virus complete genome determination and development of infectious cDNA by combining two methods- 3' RACE and RNA-Seq.

Ryegrass mottle virus (RGMoV; genus: Sobemovirus) is a single-stranded positive RNA virus with a 30 nm viral particle size. It exhibits T = 3 symmetry with 180 coat protein (CP) subunits forming a viral structure. The RGMoV genome comprises five open reading frames that encode P1, Px, a membrane-anchored 3C-like serine protease, a viral genome-linked protein, P16, an RNA-dependent RNA polymerase, and CP. The RGMoV genome size varies, ranging from 4175 nt (MW411579.1) to 4253 nt (MW411579.1) in the deposited sequences. An earlier deposited RGMoV complete genome sequence of 4212 nt length (EF091714.1) was used to develop an infectious complementary DNA (icDNA) construct for in vitro gRNA transcription from the T7 promoter. However, viral infection was not induced when the transcribed gRNA was introduced into oat plants, indicating the potential absence of certain sequences in either the 5' or 3' untranslated regions (UTR) or both. The complete sequence of the 3' UTR was determined through 3' end RACE, while the 5' UTR was identified using high-throughput sequencing (HTS)-RNA-Seq to resolve the potential absences. Only the icDNA vector containing the newly identified UTR sequences proved infectious, resulting in typical viral infection symptoms and subsequent propagation of progeny viruses, exhibiting the ability to cause repeated infections in oat plants after at least one passage. The successful generation of icDNA highlighted the synergistic potential of utilizing both methods when a single approach failed. Furthermore, this study demonstrated the reliability of HTS as a method for determining the complete genome sequence of viral genomes.

Biological factors and statistical limitations prevent detection of most noncanonical proteins by mass spectrometry.

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry (MS) experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here, we leveraged recent advances in ribosome profiling and MS to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly expressed to be detected by shotgun MS at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for 4 noncanonical proteins in MS data, which were also supported by evolution and translation data. These results illustrate the power of MS to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly expressed proteins.

Integration of genomics and transcriptomics highlights the crucial role of chromosome 5 open reading frame 34 in various human malignancies.

Although some data suggest that chromosome 5 open reading frame 34 (C5orf34) plays a pivotal part in the onset and disease progression of various cancers, there is no pan-cancer investigation of C5orf34 at present. This study sought to establish the predictive importance of C5orf34 in a variety of human malignancies and to understand its fundamental immunological function. In our research, we applied a combination of several bioinformatics techniques and basic experiments to investigate the differential expression of C5orf34, and its relationship with prognosis, methylation, single nucleotide variant, clinical characteristics, microsatellite instability, tumor mutational burden, copy number variation, and immune cell infiltration of several cancers from the database that is publicly available with the aim of identifying the potential prognostic markers. In this study we found that C5orf34 expression differed significantly among cancers types, according to the findings. The expression level of C5orf34 is markedly increased in the majority of malignancies when compared to normal tissues, which is significantly correlated with an unfavorable prognosis of patients. Immunohistochemical staining confirmed the findings that C5orf34 expression was remarkably up-regulated in a variety of gynecologic cancers. Moreover, C5orf34 expression was shown to be correlated with the clinical features of patients. C5orf34 was also found to be expressed with genes that code for the major immune suppressors, chemokines, immune activators, chemokine receptors, and histocompatibility complex. Finally, our study shows that C5orf34 has the potential to be employed as a prognostic biomarker. Moreover, it might regulate the immune microenvironment in a variety of malignancies.

Genome sequence and characterization of a novel Pseudomonas putida phage, MiCath.

Pseudomonads are ubiquitous bacteria with importance in medicine, soil, agriculture, and biomanufacturing. We report a novel Pseudomonas putida phage, MiCath, which is the first known phage infecting P. putida S12, a strain increasingly used as a synthetic biology chassis. MiCath was isolated from garden soil under a tomato plant using P. putida S12 as a host and was also found to infect four other P. putida strains. MiCath has a ~ 61 kbp double-stranded DNA genome which encodes 97 predicted open reading frames (ORFs); functions could only be predicted for 48 ORFs using comparative genomics. Functions include structural phage proteins, other common phage proteins (e.g., terminase), a queuosine gene cassette, a cas4 exonuclease, and an endosialidase. Restriction digestion analysis suggests the queuosine gene cassette encodes a pathway capable of modification of guanine residues. When compared to other phage genomes, MiCath shares at most 74% nucleotide identity over 2% of the genome with any sequenced phage. Overall, MiCath is a novel phage with no close relatives, encoding many unique gene products.